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Perturbative Quantum Chromodynamics (pQCD) predicts that the small- a; gluons in a hadron 
wavefunction should form a Color Glass Condensate (CGC), which has universal properties, which 
are the same for all hadrons or nuclei. Assuming this property, in this paper we cross relate the 
l/"} ' current CGC descriptions of the ep HERA data and dAu RHIC data. In particular, we use the 

quark dipole scattering amplitude recently proposed by Kharzeev, Kovchegov and Tuchin (KKT) 
to explain the high pr particle suppression observed in dAu collisions at RHIC in our calculations 
of the proton and longitudinal structure functions. We present a detailed comparison between this 
(3JT)' parameterization and those proposed to describe the ep HERA data. We find out that, due to 

^ , its peculiar dependence on the energy and dipole separation, the KKT parameterization is able to 

■ describe the experimental ep data only in a limited kinematical range of photon virtualities. 
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In the past few years much theoretical effort has been devoted towards the understanding of the growth of the 
I/-) , total scattering cross sections with energy. These studies are mainly motivated by the violation of the unitarity (or 
Froissart) bound by the solutions of the linear perturbative DGLAP Q and BFKL Q evolution equations. Since 
these evolution equations predict that the cross section rises obeying a power law of the energy, violating the Froissart 
bound Q, new dynamical effects associated with the unitarity corrections are expected to stop its further growth 
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0,01 . This expectation can be easily understood: while for large momentum transfer k±, the BFKL equation predicts 
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that the mechanism g — > gg populates the transverse space with a large number of small size gluons per unit of 
rapidity (the transverse size of a gluon with momentum k± is proportional to l/k±), for small k± the produced gluons 



, — i 

X 



overlap and fusion processes, gg — > g, are equally important. Considering the latter process, the rise of the gluon 
distribution below a typical scale is reduced, restoring the unitarity. That typical scale is energy dependent and is 
called saturation scale Q s . The saturation momentum sets the critical transverse size for the unitarization of the 
cross sections. In other words, unitarity is restored by including non-linear corrections in the evolution equations 
SEBSESHIHEIIl- Such effects are small for k\ > Q 2 and very strong for k\ < Q 2 , leading to the 
saturation of the scattering amplitude. 

In the high energy limit, perturbative Quantum Chromodynamics (pQCD) predicts that the small- x gluons in a 
hadron wavefunction should form a Color Glass Condensate (CGC), which is described by an infinite hierarchy of the 
coupled evolution equations for the correlators of Wilson lines HHUElEl- I n the absence of correlations, the first 
equation in the Balitsky-JIMWLK hierarchy decouples and is then equivalent to the equation derived independently 
by Kovchegov within the dipole formalism [13] . A complete analytical solution of the Balitsky-Kovchegov (BK) 
equation is still lacking though there have been interesting recent developments in this direction (for recent reviews 
sec, e.g. 

0HHG3) a remarkable feature which emerges from the solution of this equation is that the dense, 
saturated system of partons to be formed in hadronic wave functions at high energy has universal properties, the 
same for all hadrons or nuclei. In particular, as the parton densities present in dAu collisions at RHIC are not 
too different from those measured in DIS at HERA, one expects CGC physics (and thus the presence of an energy 
dependent saturation scale Q s ) to affect particle production rates and cross sections. This allows us to cross relate 
these experiments in this respect and gain a clear understading of the CGC in high energy experiments. In order to 
illustrate this statement, in Fig. ^we present the A and x dependence of the saturation scale, assuming the empirical 
parameterization Q 2 = As x Qq (^t) a , with the parameters Q\ = 1.0 GeV 2 , x = 0.267 x 10 -4 and A = 0.253 as in 
Ref. [l8|. We can observe that, while in the proton case we need very small values of x to obtain large values of Q 2 , 
in the nuclear case a similar value can be obtained for values of x approximately two orders of magnitude greater. 
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FIG. 1: Saturation scale for different values of A and x. 

In particular, the value of Q 2 = 2 GeV 2 , which is estimated from ep HERA data, can be obtained in dAu collisions 
at RHIC in the forward rapidity region. A strong support for the universality of the CGC physics has been given 
recently in Ref. 0], which has noticed that the results for different collision systems in 7*p(A), dA and AA can 
be related through the geometric scaling property, which is one of the main characteristics of the high density QCD 
approaches |20L l2l| ■ 

The search of signatures for t he p arton saturation effects has been an active subject of research in the last years (for 
recent reviews see, e.g. [lj, QjJ 123 ) ■ I n particular, it has been observed that the HERA data at small x and low Q 2 
can be successfully described with the help of saturation models [1^ I24L I25I |26| . Moreover, experimental results for 
the total cross section [^3 and also for inclusive charm production [23] present the property of geometric scaling. On 
the other hand, the recently observed [2t| suppression of high pt hadron yields at forward rapidities in dAu collisions 
at RHIC has the behavior anticipated on the basis of CGC ideas [3(j ■ Although the data is qualitatively consistent 
with the predictions based on the CGC picture, only recently more quantitative analysis have been made |3lll32| (See 
also Refs. 0,133]). These approaches consider different basic assumptions in order to describe the experimental data. 
In particular, they consider distinct prescriptions for the dipole target cross section which is one of the basic elements 
of the CGC approaches. In Ref. [31| a generalization of the parameterization proposed by Iancu, Itakura and Munier 
(IIM) to describe the HERA data was used, obtaining a good agreement with the BRAHMS data on charged hadron 
production in the limited region of low transverse momenta and forward rapidity (y = 3.2). A comparison between 
this model and the RHIC data in the full kinematical range is not possible due to the behavior of the Fourier transform 
of the IIM dipole target cross section at intermediate transverse momenta |34| (For a recent detailed discussion of 
this subject see Ref. H^)- In [32|. Kharzeev, Kovchegov and Tuchin (KKT) introduced a new parameterization with 
the free parameters fitted to RHIC data. In order to describe the hadron production in dAu collisions at forward and 
mid-rapidities the authors has considered the contributions of gluon and valence quark production and convoluted it 
with the fragmentations functions and deuteron parton distributions. In particular, the gluon production cross section 
is given in terms of the gluon dipole scattering amplitude J\fc(r,x), while the valence quark production cross section 
is a function of the quark dipole scattering amplitude A/g(r,x). In principle both AfG(f,x) and J\fg(r,x) should 
be determined from the solution of the BK (or JIMWLK) evolution equation. However, as an analytical solution 
of this equation has not been accomplished so far, the authors from Ref. [23 have proposed a phenomenological 
parameterization for these two scattering amplitudes, inspired in the approximated analytical solutions of the BK 
equation for the saturation and color transparency regimes. It is important to emphasize that the pair separation 
and energy dependence proposed for the scattering amplitudes A/g(»", x) and A/q(t", x) are identical, characterized by 
a (rQ 2 ) 7 ^ Yr ' dependence, where the form of the anomalous dimension j(Y,r 2 ) is constructed considering known 
analytical solutions to the BFKL equation. As our goal in this paper is to analyze the saturation physics in the 
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deep inelastic scattering processes, which is directly associated to the quark dipole scattering amplitude, we will 
only consider the expression for Mq{r,x) proposed in Ref. j^. The main uncertainty present in this procedure is 
associated with the normalization of the dipole cross section, which comes from the impact parameter dependence, 
and is not specified in |32|. In what follows we will consider the normalization as a free parameter to be fixed in a 
comparison with the experimental data and keep all other parameters fixed as in Ref. |32| . 

Based on the universality of the hadronic wave function, we might expect that the KKT parameterization would 
also describe the HERA data on proton structure functions. The main goal of this paper is to check this expectation. 
We will compare the predictions made with the KKT cross section and HERA ep data in the kinematical region 
where the saturation effects should be present (small x and low Q 2 ). Moreover, we analyze in detail the quark dipole 
scattering amplitude proposed in Ref. [32 and compare with those previously proposed to describe HERA data 
(for a related discussion see Ref. [3(|). We observe that there are large differences in the energy and pair dipole 
size dependences of these models. We will arrive at the conclusion that the experimental data on proton structure 
function are described using the KKT prescription for the dipole cross section only in a limited kinematical range of 
photon virtualities. This result can be interpreted as an indicative that pre- asymptotic effects cannot be disregarded 
for the kinematical range of the RHIC and HERA colliders, which implies that the property of universality is still 
not manifested. It is important to emphasize that numerical studies of the BK equation show that its solution for 
intermediate rapidities presents a strong dependence in the choice for the initial condition [l3ll3^ . l3^ |. Other important 
aspect that deserves more detailed analyzes is that the impact parameter dependence, which is disregarded in the 
phcnomenological parameterizations, may have a significant effect on the behavior of the dipole scattering amplitudes. 

A comment is in order here. A systematic comparison between the IIM and KKT dipole cross sections was started 
in Refs. [3(| and 01 1 where possible ways to choose which one is the most appropriate were discussed. It was pointed 
out that one promising observable is the nuclear modification factor for photon production in deuteron-gold collisions 
at y = 3.8. In j3|| it was mentioned that the KKT parameterization had not been checked against DIS data on proton 
targets at HERA. In this paper we perform this check. 

This paper is organized as follows. In the next section we shortly review the deep inelastic scattering in the color 
dipole picture, where the relation between the proton structure function F2(x, Q 2 ) and the dipole target cross section 
&dip becomes explicit. Moreover, we review the parameterizations proposed in the literature for the quark dipole 
scattering amplitude, with particular emphasis on that proposed in Ref. [32 • In Section ITTT1 we present a detailed 
comparison between the distinct parameterizations for J\f(r,x) discussed in the previous section. The asymptotic 
predictions for the color transparency and black disk regimes are compared as well as the energy and dipole size 
dependences. Furthermore, a comparison of predictions with the F2 HERA data in the kinematical region of small 
values of x and Q 2 is presented. As a by product, we also present a comparison with the HERA data on the longitudinal 
structure function. Finally, in Section llVI we summarize our main results and conclusions. 



We start from the space-time picture of the eletron-proton/nuclei processes [3!j. The deep inelastic scattering 
ep(A) — > e + X is characterized by a large electron energy loss v (in the target rest frame) and an invariant momentum 
transfer q 2 = —Q 2 between the incoming and outgoing electron such that x — Q 2 flm^v is fixed (rnjv is the target 
mass). In terms of Fock states we then view the ep(A) scattering as follows: the electron emits a photon (|e>— > \ej>) 
with E 1 = v and p 2 ^ « Q 2 . Afterwards the photon splits into a qq (\ej>— ► \eqq>) and typically travels a distance 
l c pa l/m^x, referred to as the coherence length, before interacting in the target. For small x, the photon is converted 
into a quark pair at a large distance before the scattering. Consequently, the space-time picture of the DIS in the 
target rest frame can be viewed as the decay of the virtual photon at high energy into a quark-antiquark pair (color 
dipole), which subsequently interacts with the target (for a review see, e.g., Ref. [Hj). In the small x region, the 
color dipole crosses the target with fixed transverse distance r between the quarks. The interaction j*p(A) is further 
factorized and is given by [39j, 



where z is the longitudinal momentum fraction of the quark of flavor /. The photon wave functions ^l,t are 
determined from light cone perturbation theory and are given by 



II. DEEP INELASTIC SCATTERING 
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e 2 - {Q 2 z 2 (l-z) 2 ^ 2 (er)} 
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The variable r defines the relative transverse separation of the pair (dipole) and z (1— z) is the longitudinal momentum 
fraction of the quark (antiquark). The auxiliary variable e 2 = z(l — z) Q 2 + m 2 ? depends on the quark mass, m/. The 
K Q i are the McDonald functions and the summation is performed over the quark flavors. 

The dipole hadron cross section a dip contains all information about the target and the strong interaction physics. 
In the Color Glass Condensate (CGC) formalism [9!. IToL fTl| . a dip can be computed in the eikonal approximation and 
is given by: 



2 J d 2 bAf(x,r,b) 



(4) 



where TV is the quark dipole-target forward scattering amplitude for a given impact parameter b which encodes all 
the information about the hadronic scattering, and thus about the non-linear and quantum effects in the hadron wave 
function. The function TV can be obtained by solving an appropriate evolution equation in the rapidity Y = ln(l/x). 
The main properties of TV are: (a) for the interaction of a small dipole (r <C 1/Q S ), TV(r) ps r 2 , implying that 
this system is weakly interacting; (b) for a large dipole (r 3> 1/Q S ), the system is strongly absorbed and therefore 
TV(r) ps 1. This property is associated to the large density of saturated gluons in the hadron wave function. It is 
useful to assume that the impact parameter dependence of TV can be factorized as Af(x, r, b) = J\f(x, r)S(b), so that 
(Jdip{x,r) — &o Af(x,r), with tro being a free parameter related to the non-perturbative QCD physics. 

Several models for the dipole cross section have been used in the literature in order to fit the HERA data. Here we 
will consider only the models proposed in Refs. 0, |2j| which capture the main properties of the CGC physics. An 
equally good fit has been obtained in Ref . , where the x dependence of the dipole cross section was derived from 
the numerical solution of the BK equation, including DGLAP corrections. In Ref. j2^| Golec-Biernat and Wusthoff 
(GBW) have proposed a phenomenological saturation model where TV is given by 



Af(x,r) = 



1 — exp — 



(Qs(x)r) 



(5) 



with Q 2 = Ql e *M*o/aO. The parameters were obtained from a fit to the HERA data yielding ct = 23.03 (29.12) mb, 
A = 0.288 (0.277) and x = 3.04 • 10~ 4 (3.41 • 10 -4 ) for a 3-flavor (4-flavor) analysis [23. An additional parameter is 
the effective light quark mass, m/ = 0.14 GeV, consistent with the pion mass. It should be noticed that the quark 
mass plays the role of a regulator for the photoproduction (Q 2 — 0) cross section. The light quark mass is one of the 
non-perturbative inputs in the model. The charm quark mass is considered to be m c = 1.5 GeV. A smooth transition 
to the photoproduction limit is obtained with a modification of the Bjorken variable as, 
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Observing Eq. we notice that when Q 2 (x)r 2 <C 1, the model reduces to color transparency, whereas as one 
approaches the region Q 2 (x)r 2 pa 1, the exponential takes care of resumming many gluon exchanges, in a Glauber- 
inspired way. Intuitively, this is what happens when the proton starts to look dark. Although the GBW parameteri- 
zation gives a good description of the old HERA data, it has been ruled out by the new HERA data, with a much 
higher accuracy. This shortcoming is mainly related to the fact that this model fails to describe the Bjorken scaling 
violation and its functional form is only an approximation of the theoretical non-linear QCD approaches. 

Another CGC inpired model has been proposed to described the HERA data in Ref. [TJ|. It is based on the 
understanding of the BFKL approach in the border of the saturation region [2(|. In particular, the forward scattering 
amplitude has been calculated in both leading order (LO) and next-to-leading order (NLO) BFKL approaches in the 
geometric scaling region ^l). It reads 



TV(x,r) = [r 2 Q 2 s {x)] ls exp 



2 Q 2 S ) 



2f3a s Y 



(7) 



where the power 7^ is the (BFKL) saddle point in the vicinity of the saturation line Q 2 — Q 2 (x). In this model the 
overall normalization of the dipole cross section is given by do = 27ri? 2 , where R p is the proton radius. In addition, 
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the anomalous dimension is defined as 7 = 1 — 7 S . As usual in the BFKL formalism, a s = N c a s /n and j3 ~ 28 £(3). 
The quadratic diffusion factor in the exponential gives rise to the scaling violations, which are essential to describe 
the HERA data. As the forward scattering amplitude in Eq. (JJJ) does not include an extrapolation from the geometric 
scaling region to the saturation region, the authors from Ref. 18] have constructed a parameterization for Af(x,r) 
which smoothly interpolates between the limiting behaviors analytically under control: the solution of the BFKL 
equation for small dipole sizes, r <C l/Q s (x), and the Levin- Tuchin law f° r larger ones, r ^> 1/Q s (x). A fit 
to the structure function F2(x,Q 2 ) was performed in the kinematical range of interest, showing that it is not very 
sensitive to the details of the interpolation (for a comprehensive phenomenological analysis of the HERA results 
using the numerical solution of the BK equation see Ref. 38] ). The dipole-target forward scattering amplitude was 
parametrized as follows, 

AAo(^) KS KX¥ , for rQ s (x)< 2, (8) 
l-cxp- aln2 ( brc ^), for rQ s {x)> 2, 

where the expression for rQ s (x) > 2 (saturation region) has the correct functional form, as obtained either by 
solving the Balitsky-Kovchegov (BK) equation 0,0], or from the theory of the Color Glass Condensate (CGC) |l4j |. 
Hereafter, we label the model above by IIM. The coefficients a and b are determined from the continuity conditions of 
the dipole cross section at rQ s (x) — 2. The coefficients 7 S = 0.63 and k = 9.9 are fixed from their LO BFKL values. 
In our further calculations we shall use the parameters R p = 0.641 fm, A = 0.253, xq — 0.267 x 10~ 4 and JVq = 0.7, 
which give the best fit result. Recently, this model has also been used in phenomenological studies of vector meson 
producti on l43ll and diffractive processes 0j at HERA as well as for the description of the longitudinal structure 
function [451 ] . 

On the other hand, Kharzeev, Kovchegov and Tuchin (KKT) have proposed a new parameterization for the dipole 
scattering amplitude in order to describe hadron production in dAu collisions [32]. As already discussed in the 
Introduction, in order describe the hadron production in dAu collisions at forward and mid-rapidities these authors 
have proposed a phenomenological parameterization for the quark and gluon dipole scattering amplitudes, inspired in 
the approximated analytical solutions of the BK equation for the saturation and color transparency regimes. In this 
model the expression for the quark dipole-target forward scattering amplitude (hereafter Mq = J\f(r,x)) is given by 
HI: 



Af(r, x) — 1 — exp 
where the anomalous dimension j(Y, r 2 ) is 



-(r^Ql 

4 V N r " 



(9) 



where c is a free parameter and 



2 V £(Y>2) + v /2£(r,r2) + 7C(3)c 



ur ' rj (A/2)(y-y ) [ ' 



The authors assume that the saturation scale can be expressed by Q 2 {Y) = A 2 A 1 ^ (^) A - The form of the anomalous 
dimension is inspired by the analytical solutions to the BFKL equation 0. Namely, in the limit r — > with Y fixed 
we recover the anomalous dimension in the double logarithmic approximation 7 w 1 — -y/l/(2£). In another limit of 
large Y with r fixed, Eq. (|10fl reduces to the expression of the anomalous dimension near the saddle point in the 
leading logarithmic approximation 7 w i + 14e ^( 3 - ) ■ Therefore Eq. (|10fl mimicks the onset of the geometric scaling 
region |l8ll20l]. In the calculations of Ref. [32| it is assumed that a characteristic value of r is r w 1/(2 kx) where kr 
is the transverse momentum of the valence quark and 7 was approximated by j(Y, r 2 ) « j(Y, 1/(4 fc^)). As our goal 
is to apply this model to deep inelastic scattering, we explore two other possible approximations which are r ss l/Q s 
and r w 1/Q- In the above expressions the parameters A = 0.6 GeV and A = 0.3 are fixed by DIS data 23]. The 
initial saturation scale used in i|TT|) is defined by Q 2 — Q 2 (Yq) with Y being the lowest value of rapidity at which 
the \ow-x quantum evolution effects are essential. When applied to describe RHIC data, the amplitude Af(r,x) must 
be convoluted with the quark distribution function in the hadron and with the fragmentation function of the quark. 
Moreover the gluon contribution must be added. These procedures introduce uncertainties in the predictions, which 
can only be estimated if the formalism is applied to other processes. 
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TABLE I: Asymptotic limits of the quark dipole scattering amplitude in different models. 
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FIG. 2: Dependence of the quark dipole scattering amplitude in the squared pair separation r 2 at different values of x. 

A comment is in order here. The main goal of the IIM and KKT parameterizations is to mimick CGC physics 
in all kinematical regions. However, currently we have some theoretical control only on the asymptotic regimes 
of saturation and color transparency. Therefore we must assume some interpolation Ansatz in order to obtain a 
parameterization that may be used in practical calculations. Moreover, although both parameterizations have similar 
form, the KKT parameterization, in contrast to IIM one, includes the double logarithmic limit as well as the correct 
behavior in the saturation and geometric scaling regions. As is well known, this limit is important for large transverse 
momentum, allowing to connect the CGC physics with the DGLAP predictions. However, what is the correct linear 
limit (BFKL/DGLAP) in the kinematical regions of HERA and RHIC is still an open question. 

III. RESULTS AND DISCUSSION 

In this section we present a detailed study between the predictions of the distinct parameterizations for the quark 
dipole scattering amplitude and a comparison of its predictions with HERA data. We start presenting in Table [I] 
the asymptotic predictions for the linear regime r <C r sat and saturation regimes r 3> r sat , where r sat = l/Q s . As 
discussed in the Introduction, the critical line dividing dense and dilute regions is the saturation scale Q s , with the 
property that the smaller the x, the denser the system gets and partons start to reinteract. The basic feature of the 
GBW, IIM and KKT models is that for a given r, these models predict that the amplitudes tend to unity at small 
values of x in contrast to the linear solution which predicts a exponential growth in this kinematical region. Moreover, 
all these models predict that the system saturates early, that is for large values of x when the dipole size is larger. 
The three parameterizations present similar functional forms for the forward scattering amplitude in the two limits, 
with the IIM presenting a residual rQ s dependence in the saturation regime, but also showing saturation for large 
values of rQ s . 




7 



1 



0.8 
"*V 0.6 
0.4 



0.2 



"0 123456789 10 

rQ 

FIG. 3: Quark dipole scattering amplitude as a function of the scaling variable rQ a . 

In Fig|2l we analyze the pair separation dependence of the quark dipole scattering amplitude for different values 
of x. As expected from the previous discussion, we observe that while the GBW and IIM parameterizations present 
a similar behavior for small r 2 , the KKT one predicts a smoother dependence. In the other limit, the GBW and 
IIM parameterizations saturate for large pair separations, while the KKT one still presents a residual dependence, 
demonstrating that the asymptotic regime is only reached for very large pair separations. The characteristic feature 
which is evident in the GBW and IIM models is that the dipole cross section saturates for smaller dipoles when x 
assumes smaller values. An important aspect to be emphasized is the large difference between the predictions in 
the transition region, which we expect to be probed at HERA. For comparison we also present the predictions for 
JV from the numerical solution of the BK equation as obtained in Ref. |38(. We have that this solution has, as 
expected, the color transparency and saturation limits for small and large r, respectively. However, the transition 
region is characterized by a sharp transition in the two values of x shown. Moreover, it is important to emphasize 
the large difference between this results and the KKT prediction. In what follows we will restrict our analyzes for the 
phenomenological parameterizations which has its parameters fixed by the experimental data. 

All models have the property of geometric scaling observed in the solutions of the BK equation. Mathematically, 
geometrical scaling means that the solution of the BK equation depends only on one combined variable rQ s (x) instead 
of r and x separately, i. e., Af(r,x) = Af(rQ s (x)). In Fig. [3] we show the dependence of the quark dipole scattering 
amplitude on the scaling variable rQ s . We observe that the three dipole scattering amplitudes grow in the region of 
small values of rQ s as a power of rQ s , i.e. Af(r,x) oc (rQ s ) 27c - f/ . However, j e ff is different in each model, being 1 
for the GBW model, < 1 for the IIM model and about \ for the KKT one. This implies a different rQ s dependence 
of the dipole scattering amplitudes and dipole cross sections. Since the saturation scale drives the energy dependence 
of the dipole cross section, these models present a very distinct energy dependence. This can clearly be seen in Fig. 

where we present the x dependence of the dipole scattering amplitudes for different values of the squared pair 
separation given by r 2 = 1/Q 2 . We observe that for large Q 2 (small pair separation) the dipole scattering amplitude 
is dominated by the linear limit. Since the models have different behavior in this limit, the energy dependence is 
also different, with the GBW model presenting the strongest growth at small x. The behavior predicted by the IIM 
model is similar to the GBW one. On the other hand, the KKT model predicts the smallest growth with the energy. 
At large pair separations r > r s , which characterizes the saturation regime, the GBW and IIM models predict the 
saturation of the dipole scattering ampitude, while the KKT one still presents a growth at small values x. Basically, 
the asymptotic saturation regime is only observed for very small values of x, beyond the kincmatical range of HERA. 

The basic observable measured with a great accuracy by HERA is the proton structure function F2(x,Q 2 ) which 
is directly related with the ^y*p cross section by the following expression 

F 3 {x,<?) = -r^(a£* + *£*). (12) 
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Consequently, using the color dipole picture of DIS we can directly calculate F2 for the different parameterizations 
of the quark dipole scattering ampitudes. Similarly, we can estimate the longitudinal structure function which is 
defined by Fl(x,Q 2 ) = Q 2 /(4ir 2 a em ) x a]^ v ■ In Figs. [SJandElwe present a comparison between the predictions of 
the different models and ZEUS data ^(|. We have used that <Tdi P (x,r) — aoJ\f(x,r) with ctq as given in the GBW 
and IIM parameterizations. For the KKT parameterization we have treated <7o as a free parameter and fixed its value 
by fitting the F2 data at Q 2 — 2.7 GeV 2 . Our choice for this value of virtuality is justified by the fact that in this 
region we expect that the saturation physics should be dominant. We have tested other choices and verified that 
our main conclusion is not modified (see below). The predicitions for other values of virtualities are parameter free. 
Moreover, we have considered two different choices for the typical scale present in the process, needed to calculate 
the function 7(Y,r 2 ) in the KKT parameterization. Basically, we have assumed that r « l/Q s or r w l/<3- As we 
will demonstrate below, our predictions for F2 in the kinematical range of interest are almost identical. We consider 
only few values of the photon virtuality in the region of low Q 2 , where the saturation effects must be important. As 
expected, the GBW and IIM models describe quite well the experimental F 2 data (See Fig. |3J) . On the other hand, 
the KKT parameterization is able to describe the experimental ep data only in a limited kinematical range of photon 
virtualities around of the virtuality where the normalization is fixed. The basic aspect of this parameterization is 
that the Q 2 dependence of the proton structure function cannot be described. Furthermore, this parameterization 
predicts an energy dependence, which is smoother than observed in the data. This behavior is directly related to the 
behavior present in the dipole scattering amplitude. The curve denoted KKTq in the figure represents the results 
obtained assuming r « while in the KKT curve we assume r w l/Q s in the calculation of 7(Y,r 2 ). These two 
prescriptions differ appreciably only in the large x and/or Q 2 region. 

Finally, in Fig. [5] we present the predictions of the different models for the longitudinal structure function. For 
comparison, we also present the prediction obtained using the Altarelli-Martinelli equation and the GRV98 parameter- 
ization for the solution of the DGLAP evolution equation (for details see Ref. 45] ) . In this case we have that the KKT 
parameterization describe reasonably the few data available 01 1 similarly to the GBW and IIM parameterizations. 
However, this fact is mainly associated to the large experimental error in the current data. We believe that a future 
experimental study of the longitudinal structure function will be able to discriminate the parameterizations. 




IV. SUMMARY 



Assuming the universality of the hadron wavefunction predicted by the Color Glass Condensate formalism we can 
cross relate different experiments and gain a clear understading of the QCD dynamics at high energies. In this 
paper we have studied in detail different parameterizations of the dipole scattering ampitude proposed to describe 
the HERA and RHIC data. We have observed that these parameterizations predict distinct energy and dipole size 
dependences, mainly in the interpolation region between the linear and saturation regimes. Since the experimental 
data at HERA probe exactly this kinematical domain, a comparison with the F2 data in the region of small x and low 
Q 2 is very important, since it allows to discriminate between the parameterizations. We have concluded that the KKT 
parameterization is able to describe the experimental ep data only in a limited kinematical range of photon virtualities. 
Therefore, the scaling violations of the proton structure function, observed in the HERA data, are not reproduced 
by this model. Moreover, the KKT parameterization predicts a smoother energy dependence than that observed in 
the data. As the IIM parameterization is not able to describe the RHIC data in the full kinematical range, our result 
put in check the property of universality, present in the CGC physics, for the current kinematical range of the RHIC 
and HERA experiments. In principle, it indicate that pre-asymptotic effects as for instance those associated to the 
different initial conditions present in ep and pA collisions, cannot still be disregarded and the cross relation between 
different experiments should be made with some caution as well as the interpretation of the comparison between the 
CGC predictions and the experimental data. We believe that an unified global fit of the RHIC and HERA data could 
be useful to obtain reliable predictions for the future colliders. 
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FIG. 6: Comparison between the predictions for Fl(x,Q 2 ) of the distinct models at different values of Q 2 . 

Acknowledgments 

The authors thank M. V. T. Machado and M. A. Betemps for helpful discussions. This work was partially financed 
by the Brazilian funding agencies CNPq, FAPESP and FAPERGS. 



[1] V.N. Gribov and L.N. Lipatov, Sov. J. Nucl. Phys. 15, 438 (1972); G. Altarelli and G. Parisi, Nucl. Phys. B126, 298 

(1977); Yu.L. Dokshitzer, Sov. Phys. JETP 46, 641 (1977). 
[2] L. N. Lipatov, Sov. J. Nucl. Phys. 23, 338 (1976); E. A. Kuraev, L. N. Lipatov, V. S. Fadin, JETP 45, 1999 (1977); I. I. 

Balitskii, L. N. Lipatov, Sov. J. Nucl. Phys. 28, 822 (1978). 
[3] M. Froissart, Phys. Rev. 123, 1053 (1961); A. Martin, Phys. Rev. 129, 1432 (1963). 
[4] L.V. Gribov, E.M. Levin and M.G. Ryskin, Phys. Rep. 100, 1 (1983). 
[5] A.H. Mueller and Jian-wei Qiu, Nucl. Phys. B268, 427 (1986). 

[6] J. Bartels, Z. Phys. C 60, 471 (1993); ibid. 62, 425 (1994); J. Bartels and M. Wiisthoff, Z. Phys. C 66, 157 (1995); J. 

Bartels and C. Ewerz, JHEP 9909, 026 (1999). 
[7] A. L. Ayala Filho, M. B. Gay Ducati and E. M. Levin, Nucl. Phys. B493, 305 (1997), ibid. B511, 355 (1998). 
[8] L. McLerran and R. Venugopalan, Phys. Rev. D 49, 2233 (1994), ibid. 49, 3352 (1994), ibid. 50, 2225 (1994). 
[9] I. I. Balitsky, Nucl. Phys. B463, 99 (1996), Phys. Rev. Lett. 81, 2024 (1998), Phys. Rev. D 60, 014020 (1999), Phys. Lett. 

B 518, 235 (2001); I.I. Balitsky and A.V. Belitsky, Nucl. Phys. B629, 290 (2002). 
[10] E. Iancu, A. Leonidov and L. McLerran, Nucl. Phys. A692 (2001) 583; E. Ferreiro, E. Iancu, A. Leonidov and L. McLerran, 

Nucl. Phys. A701, 489 (2002). 

[11] J. Jalilian-Marian, A. Kovner, L. McLerran and H. Weigert, Phys. Rev. D 55, 5414 (1997); J. Jalilian-Marian, A. Kovner 
and H. Weigert, Phys. Rev. D 59, 014014 (1999), ibid. 59, 014015 (1999), ibid. 59 034007 (1999); A. Kovner, J. Guilherme 
Milhano and H. Weigert, Phys. Rev. D 62, 114005 (2000); H. Weigert, Nucl. Phys. A703, 823 (2002). 

[12] Y.V. Kovchegov, Phys. Rev. D 60, 034008 (1999), ibid. 61 074018 (2000). 

[13] M. A. Braun, Eur. Phys. J. C 16, 337 (2000); N. Armesto and M. A. Braun, Eur. Phys. J. C 20, 517 (2001). 
[14] E. Iancu and R. Venugopalan, arXiv:hep-ph/0303204 
[15] A. M. Stasto, Acta Phys. Polon. B 35, 3069 (2004). 
[16] H. Weigert, arXiv:hep-ph/0501087 

[17] J. Jalilian-Marian and Y. V. Kovchegov , |arXiv:he p-ph/0505052 
[18] E. Iancu, K. Itakura, S. Munier, Phys. Lett. B590,199 (2004). 

[19] N. Armesto, C. A. Salgado and U. A. Wiedemann, Phys. Rev. Lett. 94, 022002 (2005) 
[20] E. Iancu, K. Itakura and L. McLerran, Nucl. Phys A708, 327 (2002). 



11 



[21] J. Kwiecinski and A. M. Stasto, Phys. Rev. D 66, 014013 (2002). 

[22] V. P. Goncalves and M. V. T. Machado, Mod. Phys. Lett. 19, 2525 (2004) 

[23] K. Golec-Biernat and M. Wiisthoff, Phys. Rev. D 59, 014017 (1999), ibid. D60 114023 (1999). 

[24] J. Bartels, K. Golec-Biernat, H. Kowalski, Phys. Rev. D 66, 014001 (2002). 

[25] H. Kowalski and D. Teaney, Phys. Rev. D 68, 114005 (2003) 

[26] J. R. Forshaw and G. Shaw, JHEP 0412, 052 (2004) 

[27] A. M. Stasto, K. Golec-Biernat and J. Kwiecinski, Phys. Rev. Lett. 86, 596 (2001). 
[28] V. P. Goncalves and M. V. T. Machado, Phys. Rev. Lett. 91, 202002 (2003) 

[29] I. Arsene et al. [BRAHMS Collaboratio n], Phys. Rev. Lett. 91, 072305 (2003); Phys. Rev. Lett. 93, 242303 (2004); Phys. 

Rev. Lett. 94, 032301 (2005); arXiv:nucl-ex/0410020 
[30] R. Baier, A. Kovner and U. A. Wiedemann, Phys. Rev. D 68, 054009 (2003); J. Jalilian-Marian, Y. Nara and R. Venu- 

gopalan, Phys. Lett. B 577, 54 (2003); D. Kharzeev, Y. V. Kovchegov and K. Tuchin, Phys. Rev. D 68, 094013 (2003); 

J. L. Albacete, N. Armesto, A. Kovner, C. A. Salgado and U. A. Wiedemann, Phys. Rev. Lett. 92, 082001 (2004); E. Iancu, 

K. Itakura and D. N. Triantafyllopoulos, Nucl. Phys. A 742, 182 (2004). 
[31] J. Jalilian-Marian, Nucl. Phys. A 748, 664 (2005) 

[32] D. Kharzeev, Y.V. Kovchegov and K. Tuchin, Phys. Lett. B599, 23 (2004). 
[33] A. Dumitru, A. Hayashigaki and J. Jalilian-Marian, arXiv:hep-ph/0506308 
[34] J. Jalilian-Marian, private communication. 

[35] M. A. Betemps and M. B. Gay Ducati, Phys. Rev. D 70, 116005 (2004) 
[36] J. Jalilian-Marian, Nucl. Phys. A 753, 307 (2005) 

[37] M. A. Kimber, J. Kwiecinski and A. D. Martin, Phys. Lett. B 508, 58 (2001); E. Levin and M. Lublinsky, Nucl. Phys. 
A 696, 833 (2001); M. Lublinsky, Eur. Phys. J. C 21, 513 (2001); M. Lublinsky, E. Gotsman, E. Levin and U. Maor, 
Nucl. Phys. A 696, 851 (2001); K. Golec-Biernat, L. Motyka and A. M. Stasto, Phys. Rev. D 65, 074037 (2002); K. Golec- 
Biernat and A. M. Stasto, Nucl. Phys. B 668, 345 (2003); K. Rummukainen and H. Weigert, Nucl. Phys. A 739, 183 (2004); 
E. Gotsman, M. Kozlov, E. Levin, U. Maor and E. Naftali, Nucl. Phys. A 742, 55 (2004); J. L. Albacete, N. Armesto, 
J. G. Milhano, C. A. Salgado and U. A. Wiedemann, Phys. Rev. D 71, 014003 (2005) K. Kutak and A. M. Stasto, Eur. 
Phys. J. C 41, 343 (2005); G. Chachamis, M. Lublinsky and A. Sabio Ve ra, Nucl. Phys. A 748 649 (2005); T. Ikeda 
and L. McLerran, Nucl. Phys. A 756, 385 (2005); C. Marquet and G. Soyez, arXiv:hep-ph/0504080 R. Enberg, K. Golec- 
Biernat and S. Munier, arXiv:hep-ph/0505101 

[38] E. Gotsman, E. Levin, M. Lublinsky and U. Maor, Eur. Phys. J. C 27, 411 (2003) 

[39] N. N. Nikolaev and B. G. Zakharov, Z. Phys. C49, 607 (1991); Z. Phys. C53, 331 (1992); A. H. Mueller, Nucl. Phys. 

B415, 373 (1994); A. H. Mueller and B. Patel, Nucl. Phys. B425, 471 (1994). 
[40] V. Barone and E. Predazzi, High-Energy Particle Diffraction, Springer- Verlag, Berlin Heidelberg, (2002). 
[41] A.H. Mueller and D.N. Triantafyllopoulos, Nucl. Phys. B640, 331 (2002); D.N. Triantafyllopoulos, Nucl. Phys. B648, 

293 (2003); A.H. Mueller, Nucl. Phys. A724, 223 (2003). 
[42] E. Levin and K. Tuchin, Nucl. Phys. B 573, 833 (2000). 
[43] J. R. Forshaw, R. Sandapen and G. Shaw, Phys. Rev. D 69, 094013 (2004). 
[44] J. R. Forshaw, R. Sandapen and G. Shaw, Phys. Lett. B 594, 283 (2004). 
[45] V. P. Goncalves and M. V. T. Machado, Eur. Phys. J. C 37, 299 (2004) 
[46] S. Chekanov et al. [ZEUS Collaboration], Eur. Phys. J. C 21, 443 (2001) 

[47] C. Adloff et al. [HI Collaboration]. Eur. Phys. .1. C 21. 33 (2001); E. M. Lobodzinska, Acta Phys. Polou. B 35 (2004) 223; 
arXiv:hep-ph/0311180 



